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Abstract 

A two-nucleon potential and consistent electromagnetic currents are derived in chiral effective 
field tlieory (xEFT) at, respectively, (or N^LO) and eQ (or N^LO), where Q generically 
denotes the low-momentum scale and e is the electric charge. Dimensional regularization is used 
to renormalize the pion-loop corrections. A simple expression is derived for the magnetic dipole 
(Ml) operator associated with pion loops, consisting of two terms, one of which is determined, 
uniquely, by the isospin-dependent part of the two-pion-exchange potential. This decomposition 
is also carried out for the Ml operator arising from contact currents, in which the unique term is 
determined by the contact potential. Finally, the low-energy constants (LEC's) entering the N^LO 
potential are fixed by fits to the np S- and P-wave phase shifts up to 100 MeV lab energies. 

PACS numbers: 12.39.Fe, 13.40.-f, 21.10.Ky 
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I. INTRODUCTION, CONCLUSIONS, AND OUTLOOK 



A quantitative understanding of low-energy nuclear physics in terms of ab initio calcu- 
lations in quantum chromodynamics (QCD) is still lacking because of the non-perturbative 
character of the theory in this regime. However, the chiral symmetry exhibited by QCD 
severely restricts the form of the interactions of pions among themselves and with other 
particles In particular, the pion couples to the baryons, such as nucleons or A-isobars, 
by powers of its momentum Q, and the Lagrangian describing these interactions can be 
expanded in powers of Q/A^, where A^^ ~ 1 GeV specifies the chiral-symmetry breaking 
scale. As a consequence, classes of Lagrangians emerge, each characterized by a given power 
of Q/A^ and each involving a certain number of unknown coefficients, so called low-energy 
constants (LEC's), which are then determined by fits to experimental data (see, for example, 
the review papers [i] and [3|, and references therein). 

This approach, known as chiral effective field theory (xEFT), has been used to study two- 
and many-nucleon interactions (sj and the interaction of electroweak probes with nuclei jl, 0] . 
Its validity, though, is restricted to processes occurring at low energies. In this sense, it has a 
more limited range of applicability than meson-exchange or more phenomenological models 
of these interactions, which in fact quantitatively and successfully account for a wide variety 
of nuclear properties and reactions up to energies, in some cases, well beyond the pion 
production threshold (for a review, see Ref. [6|). However, it can be justifiably argued 
that xEFT puts nuclear physics on a more fundamental basis by providing, on the one 
hand, a direct connection between QCD and its symmetries, in particular chiral symmetry, 
and the strong and electroweak interactions in nuclei, and, on the other hand, a practical 
calculational scheme susceptible, in principle, of systematic improvement. 

Recently we derived the nuclear electromagnetic current in a xEFT with explicit pion, 
nucleon, and A- isobar degrees of freedom 0] . Formal expressions up to one loop — N^LO or 
eQ in the power counting scheme, Q generically indicating the low momentum scale, and 
e being the electric charge — were obtained in time ordered perturbation theory (TOPT) 
by employing non-relativistic Hamiltonians implied by the chiral Lagrangian formulation of 
Refs. P, [qI, 10]. An important aspect of the derivations in Ref. 0] is in the treatment of 



the reducible diagrams: recoil corrections, which arise from expanding the nucleon energy 
denominators in these diagrams, were found to partially cancel out the contributions from 
the irreducible diagrams. When applied to the nucleon-nucleon (NN) case, this approach 
removes explicit energy dependencies, and in fact leads, at least up to one loop (N^LO or 
Q^), to the same potential constructed by Epelbaum et al. by the method of the unitary 
transformation [lO|. It also generates N^LO currents, which satisfy current conservation 
with this potential. 

One-loop xEFT currents have also been derived, with nucleons and pions only, within the 
heavy-baryon (HB) formalism by Park et al. in Ref. [5|, and have been used in calculations 
of the n-p [1, 11| and n-d [12] capture cross sections, spin observables in n-p capture 13], 
and magnetic moments of the deuteron and trinucleons Tl[, by evaluating the relevant 
transition matrix elements between wave functions obtained from realistic potentials, i. e. in 
the hybrid approach. Later in the present work we shall show that there are differences 
between the currents obtained in the HB and TOPT formalisms, some of which have to do 
with the treatment of reducible diagrams mentioned above (see Sec. HV]) . We should note 



that electromagnetic currents in the isoscalar sector were also discussed in Refs. |14l . |15 
and used in calculations of the deuteron static properties and elastic form factors. 
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The N^LO (e Q°) currents, namely without loop corrections, were used in Ref. |3] to carry 
out hybrid calculations of the magnetic moments of A=2 and 3 nuclei, and thermal neutron 
radiative captures on protons and deuterons. To have an estimate of the model dependence 
due to short-range phenomena, the variation of the predictions was studied as a function 
of the cutoff parameter A, needed to regularize the two-body operators entering the matrix 
elements, as well as a function of the input potentials. These N^LO calculations did not 
provide a satisfactory description of the experimental data, particularly for the suppressed 
process ^B.{n, 7)^H, which exhibited a pronounced sensitivity to variations in A. This clearly 
pointed to the need of including loop corrections. 

This work represents the next stage in the program initiated in Ref. [3]. It constructs, 
consistently within the xEFT framework, a NN potential and one- and two-body currents 
up to N'^LO, with the ultimate aim of studying electromagnetic properties and radiative 
captures in few-nucleon systems at this order. More specifically, it fulfills two objectives. 
The first is the construction, in dimensional regularization, of a A^A^ potential at one loop 
(Sec. ini)- The nine LEC's — Cs and Ct at Q*^, and Ci,...,C7 at in the notation of 
Ref. [LQl — which enter the potential at this order are determined by fitting the np S- and 
P-wave phase shifts u p to 100 MeV lab energies, obtained in the recent partial-wave analysis 
of Gross and Stadler [l6| (Sec. |V]). Differences between the present version of the potential 
and that obtained by Epelbaum and collaborators 17| are not substantive, since they relate 
to the use of a different form for the regulator in the solution of the Lippmann-Schwinger 
equation and the adoption, in their case, of the older Nijmegen phase-shift analysis [18j] for 
the determination of the LEC's. 

The second objective is to carry out the renormalization (in dimensional regularization) 
of the tree-level and one-loop two-body currents, and to derive the complete set of contact 
currents at N'^LO (Sec. lIII|) . Those implied by minimal substitution in the contact interaction 
Hamiltonians with two gradients of the nucleon fields were in fact obtained in Ref. 0]. 
However, in that work non-minimal couplings were not considered: we remedy that omission 
here. Lastly, in the present study we also derive (renormalized) expressions for the magnetic 
dipole (Ml) operator at N^LO (Sees. [TVl and IVll) . We find it convenient to separate, in 
the contributions from loop corrections, a term dependent on the center-of-mass position 
of the two nucleons which is uniquely determined via current conservation by 

the isospin-dependent part of the two-pion-exchange chiral potential, and a translationally- 
invariant term. The latter has a different isospin structure than that of Ref. [sj for the 
reason mentioned earlier. 

This decomposition is carried out also for the Ml operator generated by the N^LO contact 
currents. The center-of-mass dependent term is related to the contact potential, specifically 
the part of it which is momentum-dependent and therefore does not commute with the 
charge operator. However, the translationally invariant contact Ml operator depends on 
two LEC's. There are also N^LO (translationally invariant) Ml corrections at tree level, 
involving one-pion exchange, which depend on three additional LEC's. These five LEC's 
could be fixed either by reproducing a combination of nucleon and nuclear data — for example, 
pion-photoproduction data on a single nucleon along with the deuteron magnetic moment 
and cross section for np radiative capture at thermal energies — or by relying exclusively on 
nuclear data — by fitting, in addition to the observables mentioned earlier, also the trinucleon 
magnetic moments and radii. In this respect, we note that there appear to be no three-body 
currents entering at N^LO (namely, eQ~^ in A=3 systems) \2l\. 

The stage is now set for carrying out a consistent xEFT calculation of electromagnetic 
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properties and reactions in A=2-4 nuclei. The thermal neutron n-d and n-^He and keV p-d 
captures are especially interesting, since the Ml transitions connecting the continuum states 
to the hydrogen- and helium-isotope bound states are inhibited at the one-body (LO) level. 
As a result, the cross sections for these processes are sign ificantly enhanced by many-body 
components in the electromagnetic current operator [l^, 2^. Work along these lines is in 



progress. However, it remains to be seen whether the N^LO operators derived in this study 
will reduce the sensitivity to short-range physics found in the N^LO hybrid calculations 
(for the n-d case) referred to earlier, and bring theory into satisfactory agreement with 
experiment. 



II. NN POTENTIAL AT N^LO 

This section deals with the construction of the A^A'^ potential in xEFT up to order Q^, or 
N^LO. It is derived by retaining only pions and nucleons as degrees of freedom — the inclusion 
of explicit A-isobar degrees of freedom is deferred to a later work j2l|. The formalism as 
well as the techniques we adopt have already been described in Ref. 0], and we will not 
reformulate them here. 

In Fig. [T]we show the diagrams illustrating the contributions occurring up to N^LO. At 
LO (Q^) there is a contact interaction, panel a), along with the one-pion-exchange (OPE) 
contribution, panel b). At N^LO we distinguish among three different categories, which 
are: i) contact interactions involving two gradients acting on the nucleons' fields, panel c); 
ii) two-pion-exchange loop contributions, panels d)-f); and iii) loop corrections to the LO 
contact interaction, panels g) and i), and to the OPE contribution, panel h). Note that in 
the figure we display only one among the possible time orderings. 



p/ -p/ 




c) d) e) /) g) h) i) 



FIG. 1: Diagrams illustrating contributions to the A^A^ potential entering at LO {Q^), panels 
a) and b), and N^LO (Q^), panels c)-i). Nucleons and pions are denoted by solid and dashed 
lines, respectively. The filled circle in panel c) represents the vertex from contact Hamiltonians 
containing two gradients of the nucleons' fields. Only one among the possible time orderings is 
shown for each contribution with more than one vertex. 

The time ordered diagrams of panels a)-e) are irreducible, while those of panels f )-g) have 
both reducible and irreducible character. In order to avoid double counting of the reducible 
contributions due to insertion of the LO potential into the Lippmann-Schwinger equation, 
the NN potential is defined as the sum of the irreducible diagrams only. 
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The evaluation of the A^A^ potential is carried out in the static limit. Corrections to 
this approximation arise from kinetic energies of nucleons, and are referred to as recoil 
corrections. The latter are not accounted for in the Lippmann-Schwinger equation in which 
only the static potential is iterated. Hence they have been included below along with the 
irreducible contributions. The resulting potential is in agreement with that obtained by 
Epelbaum et al. in Ref. [13] • Special treatment is reserved for the diagrams of panels h) 
and i), which are discussed later in this section. 



A. Formal expressions 

In what follows we use the notation introduced in Ref. 0]. In particular, the potential is 
obtained in the center-of-mass frame where the nucleons' initial and final relative momenta 



are p and p', respectively. We also define k = p' — p, K = (p' + p)/2, Uk = yk"^ + m^, and 

In the remainder of this section we will refer to the panels in Fig. [H 

The diagram illustrated by panel a) gives rise to the LO order contact potential v*^"^^, 
which is expressed in terms of two LEC's Cs and Ct as 

t;CTO = Cs + Ct(T,- (72 , (2.2) 

while that of panel b) leads to the standard OPE potential, 

^ C^) = ■ —2 • (2.3) 

Next we consider the contributions arising from panel c). There is a number of contact 
Hamiltonians involving two gradients acting on the nucleons' fields compatible with the 
required symmetries of the underlying theory. In fact, the list of fourteen given in Ref. l3 



and reported in Appendix D of Ref. |7|] is redundant, since relations exist among the terms 
proportional to C4, Cg, and Cg, and those proportional to C7, Cg, C[q, C[i (see Appendix 
We will not enforce these in the following, since, in any case, the contact Hamiltonians (all 
twelve of them) lead (in the center-of-mass frame) to seven independent operator structures 
in the potential, each multiplied by a coefficient which is a linear combination of LEC's. 
Specifically, 

v^'^'ik.K) = CiA;=^ + C2ir' + (C3A:2 + C4ir')^Ti-o-2 + zC5^^-^-Kxk 

+ Ce • k 0-2 ■ k + o-i ■ K 0-2 ■ K , (2.4) 
where the Cj's (i = 1, . . . , 7) are linear combinations of the Cj"s {i = 1, . . . , 14), given by 

Ci = C[ — C'^ + C2/2 , 
C2 = 4C[-4C'^-2C'2 , 

C3 = C'g + C'12/2 — , 

C4 = 4C^-2Cl2 + 4C;4 , (2.5) 

C, = 2C',-AC',-2C', , 

Cq = Cj + Cg + C[q/2 + C'ii/2 — Cjg , 

C7 = 4a, + 4C's-2C[o-2C[, + 4C[s , 
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and, as per the comment above, one should keep in mind that any single one of the terms 
proportional to C4, C5, and Cg (C7, Cg, C[q, and C[i) can be reduced to a combination of 
the remaining ones by a simple redefinition of the LEC's. 

The two-pion-exchange loop diagrams of panels d)-f ) generate the following contribution: 

V'^ik) =V^{k)+V,{k)+Vf{k) 



where uj± = y(p ± k)^ + Am\. Note that recoil corrections to the reducible box diagrams 
have been included in the expressions above (for a detailed discussion of this aspect of the 
present study, see Section VI of Ref. [7|]). These recoil terms need also be accounted for 
when dealing with loop corrections to the LO contact and OPE interactions. The resulting 
contributions, panels g) and h), are then found to be 

- - ^^^^CTT,-T2(T,-aJ 4, (2.7) 



= -^^--k.-k/ (2.8) 

The potential constructed so far is in agreement with that obtained by Epelbaum et al. in 
Ref. |10| by the method of unitary transformations, but for an overall factor of +8/3 rather 
than -1/3 in Eq. fl^ . 

Lastly, we consider the diagram illustrated in panel i), which has both reducible and 
irreducible parts. The former describe interactions involving "dressed nucleons" . We do not 
take into account recoil corrections arising from the pion emitted and reabsorbed by the 
same nucleon. Hence, for diagram i) we retain the irreducible part only, and obtain 

2 2 
V, = ^(3Cs-CTa,-a2) f ^ . (2.9) 

^TT "'P 

Again, this approach leads to a result which differs from that reported in Ref. TO] for 
this diagram, specifically the term proportional to Cs in Eq. (12. 9p is absent, while that 
proportional to Ct is multiplied by —4:{g\/F^) rather than —{g\/F^). However, as it 
will become clear in the next section, these differences — for diagrams h) and i) — do not 
affect the definition of the renormalized potential, since they only lead to differences in the 
renormalization of the LEC's Cs, Ct, and qa- 



B. Renormalization 

The potential defined in the previous section contains ultraviolet divergencies which need 
to be removed by a proper renormalization procedure. In order to isolate these divergencies, 
the kernels of the N^LO contributions have been evaluated using dimensional regulariza- 
tion, and the relevant integration formulae are listed in Appendix [Bl Here we sketch the 
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regularization procedure of the various contributions, and give the final expression for the 
renormahzed AW potential. 

As an example, we discuss, in some detail, the regularization of the two-pion-exchange 
contribution of Eq. fl2.6p . In terms of the kernels L(k), I^'^'^\k) and J*-^"^(/c) defined in 
Appendix [B| it reads as 



,27r/ 



TI-T2 



8F^ 



m -8g\[l^'\k) - ei^'\k)] +4g\[j^'\k) 



^{a,xkUcT,xk),jifik) . (2.10) 



By inserting the explicit expressions of these kernels in the previous equation, we obtain 

V^'^ik) = V^'^ik) + Ti ■ T2 P2{k) + O-i ■ 0-2 - CTi ■ k 0-2 • k) Po , (2.11) 

where the renormalized (finite) part of the two-pion-exchange potential, denoted by U^'^(k), 
is given by 



^-(k) 



1 



487r2 



n-T2G{k) 



Amlil + Ag\ - 5g\) + k\l + lOgl - 23g\] 



48 

4m^ + k'^ 

(2.12) 

(2.13) 



+ -^G(A;)(A;Vi-o-2-rTi-ko-2-k 



with 



G{k) 



Aml + k'^ ^ J A ml + k^ + k 
Ami + k'^ — k 



k 



where the loop function G{k) defined here differs by a factor two from that given in Ref. 10 
The divergencies are lumped into the polynomials P2{k) (of order two) and constant Pq- 



P2{k)=- 



247r2 Pi 



4 + 22g\ - 29g\ - 9g\{2 - 5g\)(- - + 7-ln7r + ln^ 



-1^ 



l + 7g\- 9g\ - - (1 + I0g\ " 23^7^ 



3gi f 2 , , m: 

— h7-ln7r + ln 



2 , , 

---|-7-ln7r + ln^ 
e /i^ 

2 3/ ■ 



,(2.14) 



(2.15) 



where the parameter e —>■ 0^, 7 is the Euler-Mascheroni constant, and /i is the renormal- 
ization scale brought in by the dimensional regularization procedure. The isospin structure 
Ti ■ T2 multiplying the polynomial P2{k) can be reduced by Fierz rearrangement so as to 
match structures occurring in the LO v'^"^^ and N^LO v'^'^'^{k) contact contributions. Indeed, 
because of the antisymmetry of two-nucleon states. 



Tl-T2 = —2 — (Ti ■ 0-2 , 
Tl ■T2k^ = -4 (1 + CTi • (T2) - k^ 



(2.16) 
(2.17) 
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It is then seen that the terms in Pq and P2{k) renormahze Cs, Ct, Ci, C2, C4 and Cq. For 
example, the last term of Eq. (12. lip is absorbed by the redefinition, 

Ce = Ca + 5j|j.-(--+7-l.. + ln--5). (2,18) 

where the factor /i^*^ is needed because the mass dimension of the LEG Ce is c/ — 7 in space 
dimensions. Note that the renormalized Cq remains /^-independent, as becomes obvious by 
taking the logarithmic derivative with respect to /t and neglecting 0(e) terms. For ease 
of notation, we will omit the overline and tacitly imply that the LECs have been properly 
renormalized. 

The contributions in Eqs. (12. 7p . (12.90 . and (12.80 lead to further renormalization of the 
LECs Cs and Ct, as well as the axial coupling constant qa entering the LO OPE: 

v, + v, = -^Crn- T2 o-i ■ 0-2 M(3) + ^(^3Cs-Ct(t,- a^) M^'^ , (2.19) 

-.(k) = -:^n.r,^lJi^iJiAfP), (2.20) 

where the constants M^"'^ are listed in Appendix [Bl The complete NN potential up to N^LO 
included is then given as 

v{k, K) = IJC™ + I;"(k) + v^^^ik, K) + I;2"(k) , (2.21) 

where w^™, v"", if^^, and v"^^ are defined in Eqs. ((221), (ESD, (EJD, and (l212ll . respectively, 
and the overline indicates that the LECs qa and some of the C'^ have been renormalized. 



III. ELECTROMAGNETIC CURRENTS 

In this section we construct the electromagnetic current operator for a two-nucleon system 
in xEFT. In the power-counting scheme of Ref. 0], the LO term results from the coupling 
of the external photon field to the individual nucleons, and is counted as eQ~^, where a 
factor e Q is from the jNN vertex, and a factor follows from the momentum 5-function 
implicit in this type of disconnected diagrams, see panel a) of Fig. [21 Higher order terms 
are suppressed by additional powers of Q, and formal expressions up to N^LO, i.e. eQ, have 
been derived in Ref. [7| . In this section, we proceed to regularize the loop integrals entering 
these N^LO currents, and to derive the corresponding finite parts. 

At this order, we distinguish among four classes of contributions: i) currents generated 
by minimal substitution in the four-nucleon contact interactions involving two gradients 
of the nucleons' fields, as well as by non-minimal couplings (these were not considered in 
Ref. 0]); ii) two-pion exchange currents at one loop; iii) one- loop corrections to tree- level 
currents; and iv) (Q/M)^ relativistic corrections to the NLO currents resulting from the 
non-relativistic reduction of the vertices. The latter are neglected in the present work. 

For completeness, we report below the expressions up to N^LO derived in Ref. 0], and 
shown in Fig. [2l As emphasized earlier, in the present study we do not explicitly include 
A-isobar degrees of freedom. In the following, the momenta are defined as 

k, = p[-pi, K, = (p: + pO/2, (3.1) 
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Pi + ki Pi + ki P2 + k. 




Pi P) Pi P2 



a) b) c) d) 

FIG. 2: Diagrams illustrating one- and two-body currents entering at LO (eQ~^), panel a), NLO 
{eQ~^), panels b) and c), and N^LO (eQ^), panel d). Nucleons, pions, and photons are denoted 
by solid, dashed, and wavy lines, respectively. The square represents the (Q/mAr)^ relativistic 
correction to the LO one-body current. Only one among the possible time orderings is shown for 
the NLO currents. 



where pi and are the initial and final momenta of nucleon i. 
The LO contribution of panel a) in Fig. [2] is 



jLO 
Ja 



2mj\[ 



+ 1 



2eN,i Ki + ifiN,i CTi X q 

where q is the photon momentum, q = kj, and 

Cat = (1 + r^)/2 , kn = {i^s + i^v'rz)/2 , = cn + hn 



(3.2) 



(3.3) 



Ks and Ky being the isoscalar and isovector combinations of the anomalous magnetic mo- 
ments of the proton and neutron (^5 = —0.12 fi^ and Ky = 3.706 ^in). Loop corrections 
to the one-body current above, occurring at NLO and N^LO, are absorbed into ns and Ky 
The NLO seagull and pion- in- flight contributions, represented in panels b) and c), are: 



•NLO 

Jb 



•NLO 
Jc 



-«e-^ (Ti X T2)^ CTi 



(T2 ■ k2 



LU 



2, 



k2 



Ze— (n X T2 



„ o-CTi ■ ki (72 ■ k2 

^4 ^4 



(3.4) 
(3.5) 



where the momenta transferred to nucleons 1 and 2 add up to q, ki + k2 = q. Lastly, the 
N^LO (relativistic) correction to the LO current, represented in panel d), reads: 



•N^LO 

Jd 



1 m 



-Cat,! 



N 

i e 

o — r^^.i 



2 (^Kl + gV4) (2 Ki + z 0-1 X q) + Ki ■ q (q + 2 z o-i X Ki) 
Ki ■ q (4 CTi X Ki - i q) - (2 i Ki - CTi X q) 



+2 (Ki X q) cTi ■ Ki 



+ 1^2 . 



(3.6) 



In addition to the classes mentioned earlier, there are N'^LO contributions [2J] involving 
the standard ttNN vertex on one nucleon, and jnNN vertices of order eQ^ on the other 
nucleon, derived from the following interaction Hamiltonian 25|] 



H. 



— / dx Art(x) \d^ V7r,(x) + rf^ Ta V7r,(x) 



(2) 

'ynNN rp 

- d'^i ezab TaO- X V7r5(x)l A^(x) ■ V X A(x) 



(3.7) 
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where the notation and conventions of Ref. Q have been adopted for the nucleon (A^), pion 
fvTa), and photon (A) fields, and d^, dg, and are related to the original couplings given 



by Fettes et al. [25| via = 8 [dg + 5fA/(64m^)] and similarly for dg, and (igi = 2(^21 + d: 



22- 



The resulting two-body current is given by 



•N^LO _ • 9A_ 
Jtree ~ ^ ^ p2 



X q ■ 



0-2 ■ k2 



a; 



+ 1 



(3.8) 



fc2 



and in principle the unknown EEC's rfg, rfg, and dl^i could be determined by pion photo- 
production data on a single nucleon or nuclear data (as discussed in Sec. H]). The isovector 
part of jj^ee ^ the same structure as the current involving iV- A excitation 0] , to which it 
reduces if the following identifications are made: d^i/dg = 1/4, and d^ = 4/x*/i^/(9 m^v A), 
where Ha is the vrA^A coupling constant, /i* is the A^A-transition magnetic moment, and A 
is the A-A^ mass difference, A = m/^ — rriN. 

Configuration-space representations of the current operators follow from 



f f e*k,.{r;+rO/2g^Ki.(ri-r,)^(ki-q) jW(kl,Ki) 



j(2)(q) = f [ ^^^-^^^^--^^ + 

,iq.R [ gik.rj(2)(q^]^) ^ 



(3.9) 



(3.10) 



where j*^^-* or j*^^) denote any one-body or two-body operators, respectively, and 5{. . .) = 
(27r)^5(. . .). Note that — iV^5(r^ — rj), i.e. the configuration-space representation of 

the momentum operator, and in the second line of Eq. (13.101) R and r are the center-of-mass 
and relative positions of the two nucleons. 



A. N'^LO currents: terms from four-nucleon contact interactions 

The N^LO currents obtained by minimal substitution in the contact interactions involving 
two gradients of the nucleons' fields have been constructed in Ref. [3| , and are reported below 
for reference: 



■N^LO 
JCT7 



-eei 2 {2C[-a^) ^2 + ^C'.,\^i + iC'^ (0-1 + 0-2) x k2 + ?C^o-i x ki 

- z 0-2 X ki + 2 (2 - CJo) (K2 ■ 0-2) o-i + 2 (2 - C[^) (K2 ■ 0-1) (t^ 

- 2 [(Ki ■ <Ti) (72 + (Ki • (72) o-i] + 2 (2 - C[^) K2 (0-1 ■ 0-2) 

- 4C(4Ki (0-1 -era)] +1^2 . (3.11) 



In addition to these, there are contributions due to non-minimal couplings, as derived in 
Appendix \^ 



.N3LO _ _ • 

JCT7nm — * ^ 



Cl5 CTi + C(g (ri,2 - T2,z) CTl 



X q+ 1 



(3.12) 
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Pi + ki p2 + k2 



• 1 



Pi , P2 

a) 



2 

1 




^^2 




2 . 
y 


/ 

/ 

/ 




/ 

1 


> ' 3 

? 








c) 


2 

3 


1 








5) 



h) 



FIG. 3: Diagrams illustrating one-loop two-body currents entering at N^LO (eQ), notation as in 
Fig. [2j Only one among the possible time orderings is shown for each contribution. 



B. N^LO currents: one-loop corrections 

Loop corrections entering at N^LO have been derived in Ref. 0], and their formal ex- 
pressions, corrected from a number of typographical errors, are listed, for reference, in 
Appendix [C] of the present paper. In Ref. j7l| , it was also shown that the one- loop currents 
satisfy the continuity equation with the two-pion-exchange potential of Sec. Ill A[ Here we 
discuss their renormalization. We start off by considering the currents (involving one and 
two pions) illustrated in panels a), d), f), g), h), and i) of Fig. [3l Those in panels b), c), 
and e) (involving three pions) are discussed in Sec. IIVI and Appendix [Dl since for them we 
only derive the magnetic dipole operators. In terms of the kernels defined in Appendix [Bl 
we obtain: 



•N^LO _ 
Ja 


. g\ 


J(0) 




TV 










jN^LO 
Jd — 


■ 9\ 

TT 




+ 






•N^LO _ 


TT 
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(3.13) 



2 ra,^ cTi X ka + (ti x Ta)^ ks 



+ 1 



(3.14) 
(3.15) 
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Ji 



N^LO 



-2^e|4n,, J^fiq) [Cs (a, x q),- - Ct (a, x q),. 



1^2 



(3.16) 



and the currents in panels f) and h) vanish in the static hmit Insertion of the finite parts 
of the various kernels in the expressions above then gives 



vN^LO 
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tN^LO 
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87r2Fi 



G{k2) 2 T2,z CTi X k2 + (Ti X T2)z '^2 +1^2 



(3.17) 
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2 T2^z CTi X k2 + (Ti X T2) z ^2 



- 4 Ti , (T2 X k2 



+ 1^2 



I e 



I e 
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47r2 Ff 



Vi X r2)2 G'(g) cri cr2 ■ q + 1 ^ 2 , 



Ti,z G{q) (C5 CTi X q - 0-2 X q) + 1 ^ 2 



(3.18) 
(3.19) 
(3.20) 



and the loop function G is defined in Eq. f l2.13l) . The divergent parts of the kernels lead to 
renormalization of some of the LEC's C[. They are given by 
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Att^F^ 



- 2 T2, CTi X k2 - (Ti X T2)z ks 



2 T2^z cTi X (k2 - q) - -(ti X r2)^ k2 



+ . . . ) (ti X T2)z (ki - k2) 



1^2, (3.21) 
1^2, (3.22) 
(3.23) 



T2,z 0-1 X (6 k2 - 4 ki) + 3(ri X T2), k2 +1^2, (3.24) 



10 T2,z CTi X ki + -(Ti X Ta)^ k2 
D 



+ 1^2 



(Ti X T2)^ C7 



(T2 (Ti ■ q - fTi fT2 ■ q 



Ct cr2 X q — CTi X q] + 1 ^ 2 



(3.25) 
(3.26) 
(3.27) 



where the dots denote finite contributions depending on the renormalization point. When 
combined together, all these divergencies can be absorbed by the renormalization of the 
which is not the case for the individual contributions. For instance, taking into account the 
antisymmetry properties of nucleons. 



{r2,z CTi + n^z 0-2) X q = -{ti^z cti + 1"2,^ cr2) x q 

1 



X '^2) 



(Ti (72 ■ q - (T2 (Ti ■ q 



-2^Ti,z - T2,z) (cTi - 0-2) X q 



(3.28) 
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leading to renormalization of C[q, and 

(n X T2), (k2 - ki) = -2 2 ei (1 + CTi ■ o-2)(Ki - K2) + 1 

leading to renormalization of C3, (7(4, (2C( — C'2) and (2C9 — C[2). 



2 , 



(3.29) 



C. N^LO currents: one-loop corrections to tree-level currents 

Contributions in this class are illustrated by the diagrams in Figs.lHandO After including 
all possible time orderings, we find for those in Fig. HJ 
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FIG. 4: Diagrams illustrating loop corrections to tree-level two-body currents, notation as in Fig. [2l 
Only one among the possible time orderings is shown for each contribution. 
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type b) 
type c) 
type d) 
type e) 
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•NLO 
Jb 



(3.30) 

2 , (3.31) 
(3.32) 
(3.33) 
(3.34) 
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type f ) 



type j) 
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(T2 ■ k2 



+ 1 



UJ 



k2 



Jlf{q) (k2 X q)j + l 
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type k) = type 1) 
while for those in Fig. [S] 

type m) + type n 

type o) + type p 
type q 
type r 

type u) + type v 
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TT 

•NLO 
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3F2 



(3.35) 
(3.36) 
(3.37) 

(3.38) 
(3.39) 
(3.40) 
(3.41) 
(3.42) 



where jb^° and j^^*"* are the seagull and pion- in- flight currents of Eqs. (13.41) and (13.51) . 
and the constants M'^"\ and kernels Iif{q) and Jij\q) are given in Appendix [Bl The 
contributions associated with diagrams of type h), i), s), and t) vanish, since the integrand 
is an odd function of the loop momentum p. Lastly, diagrams of type g) are of order e 0], 
and therefore beyond the order under consideration in the present study (eQ), and only a 
subset of the irreducible diagrams is retained in the evaluation of the type j) contribution, 
see Appendix [El 

A few comments are now in order. Firstly, the evaluation of the diagrams in the last 
row of Figs, m and is carried out by including recoil corrections to the reducible diagrams 
of corresponding topology. Cancellations occur between the irreducible and these recoil- 
corrected reducible contributions. This is discussed in Appendix [El 

Secondly, diagrams like those shown in Fig. [6] have not been considered since they are, 
like diagram g) in Fig. [H, of order e Q^, as can be easily surmised by using the counting rules 
given in Ref. [Tf. 

Thirdly, the contributions of type a), c), e), k)-l), m)-n), and u)-v) lead to (further) 
renormalization of qa, while those of type b) and o)-p) renormalize the pion mass, namely 
= m^(l -|- M^^^ /F!^). Thus, both types are accounted for in the (renormalized) seag- 
ull and pion- in- flight currents. Diagrams j) and r) generate form-factor corrections — their 

(2) (2) 

finite parts follow from the I-j and J^j kernels — to the nucleon and pion electromagnetic 
couplings. However, the contributions of diagrams d) and f) exactly cancel out. 



IV. MAGNETIC MOMENTS FROM PION EXCHANGES AT N^LO 

To begin with, it is worthwhile making some general considerations. The magnetic mo- 
ment operator fx due to a two-body current density J(x) can be separated into a term 
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m) n) o) 




u] v) 



FIG. 5: Same as in Fig. [H 



FIG. 6: Diagrams illustrating N'^LO (eQ^) 
the present study, notation as in Fig. [2l 



loop corrections to tree-level currents not included in 



de pen dent on the center-of-mass position R of the two particles and one independent of 
it [19[, as 

/i(R,r) = ^ R X y"dx J(x) + y"dx(x- R) X J(x) , (4.1) 

where, because of translational invariance, J(x) is actually a function of J(x — R, r), r being 
the relative position of the two particles, see Eq. fl3.10p . The first term in square brackets 
can be related via the continuity equation to the commutator of the charge density operator 
with the two-nucleon potential |1^], assumed to be of the form ti ■ T2 V"(r) but otherwise 
velocity independent (for example, the one- and two-pion-exchange potentials derived in 
Sec. ini), while the second term can be written in terms of the Fourier transform of J(x), 
denoted by j(q). We find: 



/x(R,r) 



e (ti X Ta)^ V^(r) R x r + i Vg x j(q) |q=o 



(4.2) 
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which, for our purposes, is more conveniently written in momentum space as 



/x(R,k) 



e (ti X T2)z R X Vfc v{k) + Vg X j(q, k) |q=o 



(4.3) 



where t ^fk) denotes the Fourier transform of V^(r). The first term above is Sachs' contri- 
bution |20|, denoted as /xsachs, to the magnetic moment: it is uniquely determined by the 
potential between the two nucleons. 

Therefore, the currents a)-e) in Fig. [3] generate a Sachs' magnetic moment — currents g) 
and i) do not contribute to it — given by 



7ia?(R, k) = -- e (n X R X V. (fc) , 

where vl'^{k) is the term proportional to Ti ■ T2 in Eq. (12.121) . i.e. 
1 



487r2 Fi 



G{k) 



48 

Ami + k'^ 



4m^(l + Agi - hg\) + k\l + IQg^ - 23g\) - 



(4.4) 



(4.5) 



The relation (14.41) can easily be verified by direct evaluation of (R/2) x ja-c(q = 0,k). 
The currents a)-e) and i) also generate a translationally invariant contribution, namely the 
second term in Eq. (14.31) . which reads (see Appendix [D] for details) 



yj-N^LO 



(k) 



where the functions Fi{k) are 



Foik) ai-F2{k) 



kcTi • k 

A;2 



2 

2 TT'^J:'^ 



Fo{k) = l-2g\ + 



89a 
/c^ + 4 ml 



+ G{k) 



F,{k) = 2-Qgl + p^ + G{k) 
k'^ + A mi 



4(l + c^2)^2 lQg\ml 



+ A 

A{l + ^g\)ml 



+ 



+ 



(P j^4rniy 
16 gl ml 



fc2 + 4m2 {k^ + Amiy 



(4.6) 

(4.7) 
(4.8) 



It is interesting to note that the constant 2 — 65'^ in F2{k) would lead to a long-range 
contribution of the type [t2^z (""i ■ V)V + 1^2] 1/r in the magnetic moment, which is, 
however, fictitious in the present context of an effective field theory valid at low momenta — 
in performing the Fourier transform, the high momentum components are suppressed by the 
cutoff Ci^ik). 

We now compare the magnetic moment operator derived here with that obtained in 
Ref. @. Firstly, we note that the Sachs term is ignored in that work. Of course, it vanishes 
in two-body systems because of its dependence on R. However, in A > 2 systems the center- 
of-mass position of a nucleon pair will not generally coincide with that of the nucleus, and 
therefore this term will contribute. 

Secondly, the treatment of the box diagrams, panels d) and e) in Fig. [31 is different in our 
approach, since the expressions listed in Eqs. (1D2I) and (1D7I) result from combining recoil- 
corrected reducible and irreducible diagrams. In particular, had we retained only the latter, 
the isospin structure of /x^^'"*~'(k) + /j,^^'"^(k) would have contained, in addition to terms 
proportional to r^^^, also the term proportional to (ti x T2)z present in Eq. (46) of Ref. j^. 

Lastly, we find that type a) and b) contributions in Fig. [3l which only consist of irreducible 
diagrams, are in agreement with the corresponding terms in Eq. (46) of Ref. j^. This is 
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easily seen by considering the Fourier transform of that equation. To this end, we first 
observe that ^ 

\l + z{z-l)p'/ml] =G{p)-2, 



dz In 



and then note that 

dre"^'''' 



dr e 



— ikr 



rcrrr 



dr 
d_ 

dr 



ipr 



[2 - Gip)] 



[2 - Gip)] 



3G{k) + kG'{k)-6 



(4.9) 



(4.10) 



[C{k)-2]a + kG'{k)^^^ , (4.11) 



where G'{k) denotes the derivative of G{k). Inserting these relations into Eq. (46) leads to 
a similar Eq. (14.61) . but with Cs and Ct taken to be zero, and 



Fo{k) ^ G{k) 1^2 - 
F2{k) ^ 2 - G{k) 



+ 4 
4: ml 



- 2 



k"^ + A ml 



The F2{k) above is the same as Eq. (14. 8 p (with qa set to zero to remove the box contribu- 
tions), while -Fo(^) differs from Eq. (14.71) by a constant, which gives rise to a zero-ranged 
operator — operators of this type were dropped in Eq. (46) anyway. 

To the magnetic moment operators of Eqs. (14.41) and (14.61) . one has to add the term of 
one-pion range originating from the current j^ee^ (Sec. IIIII) . given by 



H'tvee ~ p2 



(C?8 ^2,^ + rfg Ti ■ Ts) k - d'^^ (Ti X T2) z CTi X k 



(72 ■ k 



k^ 



mt 



(4.12) 



V. DETERMINING THE LEG'S: FITTING THE N^LO NN POTENTIAL 



We find it convenient to formulate the NN scattering- and bound-state problems in 
momentum space 26|. In the case of scattering, we solve for the fC- matrix 



Kiri{p\p) = viTi{p\p) + ^ Hdkk 

n Jo 



.2 „,JT5 



L" 



{p',k)- 



V 



k' 



(5.1) 



where /iat is the reduced mass, V denotes a principal- value integration, and the momentum- 
space matrix elements vjTl{p',p) of the potential are defined as in Eqs. (3.3) and (3.4) of 
Ref. [13], but for the factor of 27r in front of the integration over z = p' ■ p being replaced 
here by l/(87r), and the inclusion, in the present case, of an additional phase factor , 
which, for coupled channels, leads to mixing angles with signs conforming to the customary 
choice made in phase-shift analyses. 

The integral equations above are discretized, and the resulting systems of linear equations 
are then solved by direct numerical inversion. The principal-value integration is removed by 
a standard subtraction technique [23]. Once the i^-matrices in the various channels have 
been determined, the corresponding (on-shell) S'-matrices are obtained from 



S'^%p) = l + 2iiiNpK'''{p,p) ' l-2t^i^pK'''{p,p) 



-JTS 



-1 



-JTSi 



(5.2) 
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TABLE I: Values for the nucleon axial coupling constant g^, pion decay constant F^, neutral and 
charged pion masses niQ and m+, (twice) np reduced mass f^N, and he, used in the fits. 



9A 


(MeV) 


mo (MeV) 


m+ (MeV) 


2 UN (MeV) 


he (MeV-fm) 


1.29 


184.8 


134.9766 


139.5702 


938.9181 


197.32696 



from which phase shifts and, for coupled channels, mixing angles are easily determined [17 

The bound state (with JTS = 101 and L, V = 0, 2) is obtained from solutions of the 
homogeneous integral equations (2"" 



wl{p) 



POO 

/ dA;A;2^t;i°i,(p,A;)ti;^,(A;) 
Jo V 



(5.3) 



and from these, for later reference, the configuration-space 5*- and D-wave components follow 

as 

O POO 

dpp^ jL{pr)wL{p) . (5.4) 



UL[r 



2 

n Jo 



Before turning our attention to a discussion of the phase-shift fits, we note that the 
potential constructed in Sec. IIIBI needs to be (further) regularized because of its power-law 
behavior for large values of the momenta k and/or K. This is accomplished by including a 
high-momentum cutoff, which we take to be of the form 



CA{k,K)=e 



(5.5) 



so that the matrix elements of the regularized potential entering the i^'-matrix and bound- 
state equations are obtained from 



;^(k,K) = i;(k,K)CA(fc,ir) 



(5.6) 



and f(k, K) is defined as in Eq. f l2.2ip . In the following cutoff parameters A in the range 
500-700 MeV are considered. 

The LEC's Cs, Ct, and Ci, i = 1, . . . ,7, are determined by fitting the deuteron binding 
energy and S- and P-wave np phase shifts up to laboratory kinetic energies of 100 MeV, 
as obtained in the very recent (2008) analysis of Gross and Stadler 16|. The parameters 
characterizing the one- and two-pion exchange parts of the potential are listed in Tabled with 
the nucleon axial coupling constant qa determined from the Golberger-Treiman relation qa = 
QtxNnFt, / {2mN), where the vrA^A^ coupling constant is taken to have the value g^^^/^Aix) = 
13.63 ±0.20 [28!, l29|. In fact, in the OPE we include the isospin-symmetry breaking induced 
by the mass difference between charged and neutral pions, since it leads to significant effects 
in the ^Sq scattering length 3^, and therefore the OPE potential reads 



9\ 



3F2 



Tl ■ T2 



ml 



+ 



TTlj 



12 



A;2 + 



ml 



CTi ■ k cr2 ■ k 



(5.7) 

where T12 is the isotensor operator defined as T12 = Sti^zT2,z — ^"1 ■ T2, and mo and m+ are 
the neutral and charged pions masses. Finally, we note that the pion mass entering in the 
two-pion-exchange part is taken as = (mo + 2m+)/3. 

The best-fit values obtained for the LEC's are listed in Table HI] for A=500, 600, and 700 
MeV. The fitting strategy becomes obvious once the partial wave expansion of the potential 
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is carried out. In the case of spin-singlet {S = 0) channels, the contact components of the 
(partial-wave expanded) potential with JT and S = read: 

vfj'{p',p;CT0/2)^l- [\dzPj{z)\D, + D2{p''+p')-2Dsp'pz\CA{p',p,z) , (5.8) 



Stt 



where z = p' ■ p, Pj{z) 
of the LEG'S with Di 



is a Legendre polynomial, and the Di denote linear combinations 



3Ct, D2 = Ci- 3C3 - Ce + (C2 - 3C4 - C7)/4, and 



-D3 = Ci — 3 C3 — Ce — (C2 — 3 C4 — C7)/4. The cutoff function is even in z, and therefore for 

even (odd) J only and D2 {D3) contribute. In practice, Di and D2 have been determined 
by fitting the (np) singlet scattering length and effective range, and ^Sq phase shift, while 
D3 is determined by fitting the ^Pi phases. In the case of spin-triplet {S = 1) channels, the 



TABLE II: Values of the LEC's corresponding to cutoff parameters A in the range 500-700 MeV, 
obtained from fits to np phase shifts up to lab energies of 100 MeV. 





A (MeV) 




500 


600 


700 


Cs (fm2) 


-4.456420 


-4.357712 


-3.863625 


Ct (fm2) 


0.034780 


0.094149 


0.234176 


Ci (fm^) 


-0.360939 


-0.259186 


-0.268296 


C2 (fm^) 


-1.460509 


-0.934505 


-0.835226 


Ca (fm^) 


-0.349780 


-0.359547 


-0.389047 


C4 (fm^) 


-1.968636 


-1.717178 


-1.724544 


C5 (fm^) 


-0.870067 


-0.754021 


-0.695564 


Ce (fm^) 


0.326169 


0.301194 


0.348152 


C- (fm^) 


-0.727797 


-1.006459 


-0.955273 



situation is slightly more complicated. For uncoupled channels with J > 0, we write 



Ca(p',p,z) , 



.'2 



v'/j\p',p; CTO/2) = — dz Pj{z) [D, + {D, + D,) {p'^ 

- 'Pj-i{z) + Pj+i{z)] (2 Ds + D^) p'p 
while for the ^Pq channel (having JTS — Oil) 
<V(p',p; CTO/2) = - 

+ P^{z) i2Ds-Dg)p'p Ca{p',p,z) 



(5.9) 



(5.10) 



Here, the Di's denote the following LEG combinations: D4 = Cs + Ct, -D5 = Ci + C3 + 
(C2 + C4)/4, = Ce + C7/4, D, = C^ + Cs- {C2 + C4)/4, = Cg - Cv/4, and = C5. 
In terms of these, the contact components for coupled channels are given by 



/r(p',p; CTO/2) = i-£dz 



2D 



2J+1 



{p"+p')-(2Dj-Dg)p'pz 



2J+1 



+ Dg] p'p 



(5.11) 
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1 



dz 



57r J-i 



D.+ D 



2J+ 1 

Ca{p',p,z) , 



(p" + p')-(2/^7-/^9)p'p^ 



(5.12) 



CTO/2) 



1 JJ(J+1) /•! 



47r 2 J + 1 
2DsPjiz)p'p 









j\z 




Pj 





Pj^iiz)p" + Pj+^iz)p' 



Ca{p,p,z) 



(5.13) 



where L = ± is a shorthand for L = J ± 1, and the off-diagonal matrix element with 
— h is obtained from v^^{p',p; CTO/2) by exchanging p' ^ p. The parameters D4, 
and are then determined by fitting the deuteron binding energy, spin-triplet scattering 
length and effective range, and ^Si-^Di phases and mixing angle — the contributions of terms 
proportional to Dj, Dg, and Dg vanish in this channel. On the other hand, only the latter 
enter into the '^Pj=o,i.2 channels, and the associated phases can then be used to fit D^, Dg, 
and Dg. 

Results for the S- and P-wave phases used in the fits, as well as for the D-wave and 
peripheral F- and G-wave phases, and mixing angles ej=i,,..,4 are displayed in Figs. [71-[T2]up 
to 200 MeV lab kinetic energies. Effective range expansions and deuteron properties are 
listed in Table IIIII For reference, in Figs. I9HT21 following the original work by Kaiser et 
al. [31'], the phases obtained by including only the one- and two-pion-exchange [iT and v'^'^ , 
respectively) terms of the potential are also shown. These have been calculated in first order 
perturbation theory on the T-matrix, and hence are cutoff independent. 



Overall, the quality of the fits at N^LO is comparable to that reported in Refs. 17|, l37 
and, more recently, in Ref. |38|. While the cutoff dependence is relatively weak for the S-wave 
phases beyond lab energies of 100 MeV, it becomes significant for higher partial wave phases 
and for the mixing angles. In particular, the F- and G-wave phases, while small because of 
the centrifugal barrier, nevertheless display a pronounced sensitivity to short-range physics, 
although there are indications [39| that inclusion of explicit A-isobar degrees of freedom 
might reduce this sensitivity. Beyond 100 MeV, the agreement between the calculated and 
experimental phases is generally poor, and indeed in the ^Ds and ^F4 channels they have 
opposite sign. The scattering lengths are well reproduced by the fits (within ~ 1% of the 
data, see Table [In|) . however, the singlet and triplet effective ranges are both significantly 
underpredicted, by ~ 10% and ~ 5% respectively. 

The deuteron S- and D-wave radial wave functions are shown in Fig. [T3] along with those 
calculated with the Argonne Vis (AV18) potential (sOf. The D wave is particularly sensitive 
to variations in the cutoff: it is pushed in as A is increased from 500 to 700 MeV, but 
remains considerably smaller than that of the AV18 up to internucleon distances of ~ 1.5 
fm, perhaps not surprisingly, since this realistic potential has a strong tensor component at 
short range. The static properties, i.e. D- to S-state ratio, mean-square-root matter radius, 
and magnetic moment — the binding energy is fitted — are close to the experimental values, 
and their variation with A is quite modest. The quadrupole moment is underpredicted 
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FIG. 7: (Color online) The S-wave np phase shifts, obtained with cutoff parameters A=500, 600, 
and 700 MeV, are denoted by dash (red), dot-dash (green), and solid (blue) lines, respectively. The 
filled circles represent the phase-shift analysis of Ref. 16|] . 



by ~ 4%, a pathology common, to the best of our knowledge, to all realistic potentials 
(including the AV18). 

VI. N^LO MAGNETIC MOMENT FROM CONTACT CURRENTS 

The magnetic moment due to the contact currents originating from minimal couplings 
(Sec. 1111 A|) can also be separated into a Sachs term and one independent of the center- 
of-mass position R of the two nucleons. To this end, we first note that, because of the 
gradients acting on the nucleon fields, the A^A^ contact potential contains, in addition to the 
contribution f^^^(k, K) in Eq. (12.41) . also a contribution dependent on the pair momentum 
P = Pi + P2 = p'l + P2, given by 

v^^\k, K) = i CI ■ P X k + C* (o-i • P 0-2 ■ K - CTi ■ K cr2 • P) 

+ {C; + C:(T^■(T2)P' + C;(T,■P a2-P , (6.1) 

where the C*'s consist of the following LEG combinations 

CI = C',/2 + C'J2 , 

6*2 = 2 — 2 Cg — CiQ + Cii , 

C; = -C[ + Cy2 - C, , (6.2) 

C4 = — Cg + C12/2 + , 

Q = -C',/2 - CU2 + C'JA + C[JA + C[, . 
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TABLE III: Singlet and triplet np scattering lengths (a^ and at) and effective ranges (r^ and r^), 
and deuteron binding energy {Bd), D- to S-state ratio (r/rf), root-mean-square matter radius (r^), 
magnetic moment {^d), quadrupole moment (Qd), and D-state probability (Pd), obtained with 
A=500, 600, and 700 MeV, are compared to the corresponding experimental values (a^, r^, at. and 
rt from Ref. 32\, Bd from Ref. [i^, r]d from Ref. 34], and fid from Ref. [sS^, Qd from Ref. 36]). 





A (MeV) 




500 


600 


700 


Expt 


Os (fm) 


-23.729 


-23.736 


-23.736 


-23.749(8) 


Ts (fm) 


2.528 


2.558 


2.567 


2.81(5) 


at (fm) 


5.360 


5.371 


5.376 


5.424(3) 


n (fm) 


1.665 


1.680 


1.687 


1.760(5) 


(MeV) 


2.2244 


2.2246 


2.2245 


2.224575(9) 


Vd 


0.0267 


0.0260 


0.0264 


0.0256(4) 


Td (fm) 


1.943 


1.947 


1.951 


1.9734(44) 


/"d (/"iv) 


0.860 


0.858 


0.853 


0.8574382329(92) 


Qrf (fm2) 


0.275 


0.272 


0.279 


0.2859(3) 


Pd (%) 


3.44 


3.87 


4.77 
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FIG. 9: (Color online) Same as in Fig. [3 but for D-wave phase shifts. The dash-double-dot (orange) 
line is obtained in first order perturbation theory for the T-matrix by including only the one- and 
two-pion-exchange parts of the N^LO potential. 



Incidentally, we observe that Eqs. (12. 5p and (16.21) provide a one-to-one correspondence be- 
tween the LEC's and the coefficients of the AW contact potential. 

The (conserved) current jcT7° Eq. (13.111) gives rise to a Sachs magnetic moment 



N^LCCT 
P'Sachs 



t e 

Y 



1 + 



^2„ 



Rx 



« e Ti . - r2,; 



■Rx 



(6.3) 



where the only term in Vp^"^ with a non-vanishing commutator with the relative position 



•N^LO, 
JCT7 Im 



r is that proportional to Equation (16. 3p can be easily verified by considering (R/2) x 
0). 

The Ml operator above depends on the unknown C*, which could be determined, for 
example, by fittin g A =3 bound and scattering state properties, or Ml transitions in light 
nuclei with A > 2 [40]. Instead, here we will require that they vanish, i.e. that the contact 
potential is independent of the nucleon pair momentum. To the best of our knowledge, 
this approximation has been adopted, albeit implicitly, in all studies of A > 2 nuclei based 
on xEFT potentials. In this respect, we observe that relativistic boost corrections [4l| to 
the rest-frame f*-'^^(k, K), being proportional to ~ f ^"""^ (P^/m^), are suppressed by two 
additional powers of the low momentum scale Q relative to both f '-^"'"^(k, K) and fp"'"^(k, K). 
These corrections arise from the relativistic energy-momentum relation, Lorentz contraction, 
and Thomas precession of the spins, and are of a different nature than the P-dependent 
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FIG. 10: (Color online) Same as in Fig. [U but for F-wave phase shifts. 



terms in Vp"'"^(k, K), which result from the derivative couplings in the four-nucleon contact 
Hamiltonians. 

Under the assumption above (C* = 0) and after evaluating the commutator [r , v^^^], we 
find the Sachs magnetic moment to be given in momentum space by 



MSachs (R,k,K) = 



^2„ 



■R X 



2 {C2 + C4 cTi • 0-2) K - i C5 X k 



+ C7 (cTi 0-2 ■ K + cTi ■ K 0-2) 



(6.4) 



It is determined by C2, C4, C5, and C7, i.e. by the LEC's of the momentum-dependent terms 
in t>^"^^ which do not commute with the charge operator. In configuration space, K reduces 
to the relative momentum operator, and the pair correlation function S{t) is smeared over 
a length scale 1/A (A is the high- momentum cutoff introduced in Sec. |V|) . 

The R-independent contribution due to minimal couplings follows from the second term 
in Eq. (gD, 

(6.5) 



where we have used the relation Cg = — C5 implied by = 0, and have dropped a term 
proportional to (ti ^ + T2^z) (ci + (^2), since it vanishes when acting on antisymmetric two- 
nucleons states. However, the contribution due to non-minimal couplings, which only con- 
sists of translationally-invariant terms (the corresponding currents are transverse to q and 
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FIG. 11: (Color online) Same as in Fig. [H but for G-wave phase shifts. 



therefore unconstrained by the continuity equation), is given by 



N3LO,CT ^ _^ + a,)-e C[, (n,, - r2,,) (^i - a,) 



(6.6) 



Hence, the Ml operator due to minimal and non-minimal couplings is determined by two 
independent LEC's, one of which multiplies an isoscalar structure, while the other multiplies 
an isovector structure. The former (latter) could be determined by reproducing the deuteron 
magnetic moment (the cross section for np radiative capture or the isovector combination 
of the trinucleon magnetic moment). 
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FIG. 12: (Color online) Same as in Fig. [U but for the mixing angles ej. 
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APPENDIX A: N^LO CURRENTS FROM NON-MINIMAL COUPLINGS 

External currents enter into the chiral Lagrangian either by the gauging of spacetime 
derivatives (minimal coupling), or through their field strengths F^^, which transform co- 
variantly under chiral symmetry. In the case of the electromagnetic current, we have both 
isoscalar and isovector components. In the non-relativistic limit the allowed spin-space 
structures, at leading order, are 

e^jkFijN^auN N^N , (Al) 

which, by time-reversal symmetry, can only be associated with the flavor structures 1 ® 1, 
Ta ® Ta and ® 1 ±1 ® Tz), and 

Fij N^aiN N^ajN , (A2) 

which can only be associated with the antisymmetric flavor structure e^at Ta®Th. Using the 
Fierz-type identities for the Pauli matrices, 

(1)[1] = + (A3) 
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FIG. 13: (Color online) The S-wave and D-wave components of the deuteron, obtained with cutoff 
parameters A=500, 600, and 700 MeV and denoted by dash (red), dot-dash (green), and solid 
(blue) hues, respectively, are compared with those calculated from the Argonne vig potential (dash- 
double-dot black lines). 



+ {ct)[1] = (1][^) + (^][1) 



(A4) 
(A5) 



where (, ), [, ] denote the spinors (or isospinors) xL X2, xl^ X4) we are left with two operators: 



CT7 ,nm 



- / dx 



(A6) 

We also remark that the fourteen operators in the two-nucleon, two-derivative contact 
Lagrangian can be reduced to twelve, since, using partial integration, the following relation 
involving the vertices proportional to C4, C5 and Cg and to C7, Cg, CJg and C(^, can be 
shown to hold 



N^akN {V.NyVjN + N^N (V^A^) V^A^ 
{S^kS,l - 6u6,k) [N^aj, V^N N^ai V^N + {V^N^akN {V.NyaiN 



(A7) 
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APPENDIX B: DIMENSIONAL REGULARIZATION OF KERNELS 



In this appendix we report a list of general integration formulae |42l . |43j, useful to carry 
out the regularization of the various kernels occurring in the potential and current operators. 



1. Useful integrals 

We utilize the Feynman parameterization 

1 '•I 



AB 



dy 



[yA+{l-y)B]^ 



(Bl) 



and, in order to simplify the energy factors entering the kernels, we make use of the integral 
representations [44| : 



1 



TT 



d/3 



1 



2 /-"o 



dp- 



Having defined 



d'^p 



we have: 



1 



p (p2 + A)" 

p (p2 + A)- 

4 

p 



1 T{a-d/2) 
(47r)'^/2 T{a) 

1 dT{a- d/2 - 1 



^-{a-d/2-1) 



p (p2 + A)" {iTiy/' 



(47r)"'/2 2 r(a) 

1 d{d + 2) r(a-rf/2-2)^_(„_,/2_2) 



r a) 



(B2) 
(B3) 

(B4) 

(B5) 
(B6) 
(B7) 



where T{z) is the F-function satisfying zr{z) = F(z + 1), with asymptotic behavior for 
z — > given by 



T(z) 



12 



[B8) 



and 7 ~ 0.5772 is the Euler-Mascheroni constant. However, we note that, in order to 
preserve physical dimensions, a renormalization scale /j, has to be introduced, and therefore 
a factor should be understood in Eg . (|B4D . 

Finally, we use the following relations [45| to evaluate 



dx In |x — a I = X In I X — a I — 2 x + a In 



X + a 



dxx^ In |x^ 



(dxx^ln |x^ — a^l 



- X In |x 
3 

1 
5 



X — a 

2 



- x^ - 2 a^x + In 
3 



X + a 



X — a 



2 2 

^ ( x^ In |x^ - a^l x^ a^x^ - 2 a^x + In 



X + a 



X — a 



(B9) 
(BIO) 
■ (Bll) 
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2. Regularization of the kernels 

As an example, we sketch the regularization of the kernel I^^\k), given by 



/W(A;) 



d/3- 



UJ+UJ_{UJ++UJ_) n Jp Jo (CJ^ +/?2)(cj2 +^2) 



(B12) 



where uj± = y (p ± k)^ + 4m2. Using the Feynman integral parameterization of Eq. ( IBip 
with A = 00! + (3"^ and B = + f]"^, we obtain 



/W(A;) 



TT Jp JO 
1 



dy 



2 Jp Jo 



df3 [[p + (2 y - 1) k] V 4 [ml ~y{y~l) k'] + (3' 

-3/2 



dy p' + A [mi -y{y-l) k'] 



(B13) 



where in the second line we have also shifted the integration variable p^p+(2?/ — l)k. 
The integral over p is reduced to the form given in Eq. flB5l) with d = 3, a = 3/2, and 
A = 4 [m^ — y [y — 1) k"^]. With this choice of d and a, we are left with a F-function of 
vanishing argument. In order to isolate the divergent part of the integral, we set = 3 — e 
and study its asymptotic behavior for e ^ 0+. Using 



F 
F 



2 
3 
2 



7+ 0(e) 



TT 



An 



we find, neglecting 0(e) terms. 



J(0)(A;) 



In vr H 7 



;7r2 JO 



dy In 



yiy 



fc2 



(B14) 
(B15) 
(B16) 

(B17) 



After setting y — > (x + l)/2 and making use of Eq. ( lB9p . we obtain: 



I^'\k) 



I s s + k 2 , ,m2 ' 

— -In ^_- + ^-ln7r + ln^-2 



(B18) 



where s = ^jAm^ + k^. 
The kernels 



Pipj 



I^^\k) 
4f(^) / 

can be easily evaluated as shown above. We find: 
1 



(B19) 
(B20) 



/(2)(A;) 
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In r + 2 



k 



k 



+ 7 — In TT + In ■ 



29 



/ 2 m 
+ 18ml { — +7-ln7r + ln^ 



2 



(B21) 



k s 



s + k 



+ k' 



2 , , ml 

— +7-ln7r + ln^ 
e //^ 



+ 6mt 



+ 7 — In TT + In 



ml 5 



247r2 A;2 
Next, we note that 

from which we obtain: 



J s + A; 
A; s — k 



+ k' 



,2 



2 , , 8' 

-- + 7-ln7r + ln^ - - 



1 d 



I m„ 



. (B22) 



s + k 
k 



j f{u+,uJ)^—^ — \n — 
Jp 8n^ k s s — 



k s — k 



s ^ s + k 
k s — k 



2 , , m2 2' 

-- + 7-ln7r + ln^ - - 

^2 , , m, 

+ 7 — m TT + m 



(B23) 

(B24) 
, (B25) 



2 4^ 



1 L A;, / s , s + A; 
+ T-^4Tr^ I - In- 



1 



87r2 
5 



87r2 k"^ \k s-k 

8s^s + k , 
In r + 30 m2 



e 
2 



(B26) 



3A; 



k 



2 , , m: 
1-7 — InTT + m — ; 



29' 

45, 



- ,z . 2 , , 12^ 

+ 0^ (-T + 7-ln7r + ln-^-y 



3 V e 

The set of kernels involving the energy factor 

2 a;+ + a; 



(B27) 



2u!^uj_ {uj+ + ct;_)2 
can be reduced to those of type J*^2") (^/^j by noting that 



2^', + 



^ _ L-O ~\~ LU' 



2a;+a;_ + a;_)2 4 ip + 4 



(B28) 



and similarly for J(2)(A;), ^ (A;), j(^)(A;). 

The kernels involving the energy factor g{u!^,uj_), 

. . 3 2 u;+ + cj_ + 2u 

g[u+,u^) = : - + 



1 d 2 a;+ + a;_ 

"2 dm2 2w^ tu_ (cc;+ +u;_)2 ' 



(B29) 
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easily follow from 



1 d 



'p 8dml 
and similarly for K^'^^\k), leading to: 

1 



J(°)(A;) 



16 d(m^)2 



/W(A;) 



p 



2 , s + A; 1 
In r + 



647r2 



4 , s + A; 1 



2 + ^ 
k s s — k 



g4^2 ^2 



2 s + A; _ 1 
A; s s — A; ml 



647r2 



16 s, s + A; k'^ 



In 



A; 



+ 30 



2 m: 

1-7 — Invr + m — ; 

e /i' 



+ 



8ss + k 2 , ^ mi 8" 

8s s + A; ^k^_-^Q 
k s — k mi 



647r2 A;2 

Finally, for the kernel entering diagram e) in Fig. [H we obtain 



m 



" ^ s + A; 
6 TT^ A; s — A; 



-8ml + k 



4 2 
7 ^ + 

2 
e 



+ 7-ln7r + ln^ - - 



while for the constants M^") entering Eqs. fl^l^ - fl^:^ . 



1 



M(i) = / - 

7p tUp 8 vr^ 

^2 3 m! 



p 



2 

h7-ln47r + ln^-l 

,2 



2 , , , 

h 7 — ln47r + m 

e /i 



2 3 



(B30) 

(B31) 
(B32) 

(B33) 

(B34) 



(B35) 



(B36) 

(B37) 
(B38) 



APPENDIX C: ONE-LOOP TWO-BODY CURRENTS 



In this appendix we list the expressions for the one-loop currents derived in Ref. 
Referring to Fig. [3|, we have: 



type a) 



-2i 



egA f 2r2,^ (cTi X qg) 

7r 



[Tl X T2)^q2 



a;ia;2(a;i + 102) 



+ 1 



(CI) 
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type b) = 2i- 



egx /"qi-qs 



UJ1+UJ2 + UJ3 



J UJi UJ2 UJ3 {UJI + UJ2) {UJI + UJs) {UJ2 + ws) 



type c 
type d 

type e 

type g 
type i 



2 T2^z CTi ■ (qi X q2^ 



+ 1 



2Fi 



qi - qa uj2{uji + uj2 + ujs) - 3 Ui 



-2i 



eg\ f UJI + UJ2+ 0J1UJ2 



J ujIujI{uji+U2) 



LUi Us {uJi + Co'2) (c^l + Us) {U2 + Us) 

Ti X T2)z q2 (qi • q2) + 2 T2,z qi ■ q2 (CTI X q2 



(C2) 

(C3) 



+ 2 ri_^ q2 (72 ■ (qi x qa) 



+ 1 



2i 



4 



(qi - q3)/(^i,^2,t^3) (ti X Ta)^ (qi ■ q2)(q2 ■ qs) 



+ 2 T2,z iq.2 ■ qs) CTi ■ (q2 X qi) + 2 Ti,^ (qi • q2) (72 ■ (qs X q2) 



2i 



F3 



T 



Ti X T2 



qi - q2 cj^ + cui U2 + U2^ 



qi - q2 + Ui U2 + u\ 



u\ ul 



Ct(T2 ■ (qi X q2) 



Ul + U2 

+ 1^2, 



Ul + ^012 

Cs CTi • (qi X q2 



CTi ■ q2)(cr2 ■ qi) 



(C4) 

(C5) 
(C6) 



(C7) 



where the q^'s and Ui = {qf + m^)^/^ denote the momenta (with the flow as indicated in the 
figure) and energies of the exchanged pions, and the integration is on any one of the q^'s, 
the remaining qj's with j ^ i being fixed by momentum-conserving 5-functions. Lastly, the 
function f{ui,U2,Us) in the type e) current is defined as 



f(ui,U2,Us) 



+ 



Ul U2 Us{Ui + U2){Ui + Us){u2 + Us) 
{Ui + U2) {U2 + Us) {uf + Us) U2 



uf UI2 Us 



+ 



Ul U2 + U2US + Ul Us 
Ul U2 Us 

Ul Us ul 



(CJ 



APPENDIX D: MAGNETIC MOMENTS FROM LOOP CURRENTS 

In this appendix we list the translationally invariant contributions to the magnetic 
moment — second term in Eq. (14.31) — associated with currents a)-e) and i) in Fig. [31 The 
contributions of currents c) and g) vanish, while those of currents a), d), and i) read: 



7ir^°(k) 



7ir^o(k) 



+ 



+ 



8 7r2F4 
4 



r2Mk) 



2 mi 



^2„ 



o-i 



\ Ami. + k"^ 
k cTi ■ k 

A;2 



2 mi 



k cTi ■ k 



4 m?. + k"^ k"^ 



87r2 

2 ml 



r2,zG{k) 




Aml + k^ (4m2+/c2)2 



+ 



Aml + k'^ {Aml + k"^)' 



kcTi ■ k 

k^ 



egA 

87r2 Fi 



^2,2 



Ami 



Ami + k"^ 



[Dl) 
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- 1 



4 mz 



k cr^ • k 



4 ml + k'^ J kP' 



1 ^ 2 



27r2 F2 



ri,z {Cs(Ti-Ct(T2) + 1^2 



(D2) 
(D3) 



Finally, in terms of the kernels J^"^ and K'^'^\ the contributions resulting from currents 
b) and e) are given by 



2Fi 



4'^ (A;) - h kj J^''\k)] (Ti,, - [j(2)(A;) - /tV^^H'^) 



+ 1^2 



(D4) 



Mr^°(k) 



^2,2 



- k^Kg\k) - kikjK^^\k) + kik.jeK^^\k)\ en,,- 



+ 1^2, (D5) 



from which the renormalized operators follow as 

2m! \ 
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2 k CTl • k 

4m2 + A;2 /^^ 



e k CTl ■ k 



ivr- 



T2 
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APPENDIX E: RECOIL CORRECTIONS 

Consider the set of time-ordered diagrams, displayed in Fig. [H] and denoted as type i) in 
Fig. m It is easily seen that recoil corrections in diagrams a)+b) and i)-l-j) cancel out the 
contributions associated with diagrams c)+d) and k)+l), respectively, so that the expression 
for type i) diagrams in Fig. 0] — which happens to vanish — results from diagrams e)-h). Let N 
denote the product of the four vertices in diagrams a)-d); then the contribution of diagrams 
a)+b) is given by 



a) + b) of Fig. 14 



N 



{Ei -E'-E2 + tv){Ei ~Ep-E2-iUi + ir]) 



X 



Ei- E[- E2-UJ2 + ir] Ei- E'- E'2- UJ2 + ir] 



, (El) 
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FIG. 14: Set of time-ordered diagrams for the contribution illustrated by the single diagram i) in 
Fig. m Notation as in Fig. El 



where the labeling of the momenta is as in panel a), and Ep and E'^ are the energies of the 
intermediate nucleons. The expression in square brackets above can be expanded as 





1 






U!2 



2 + 



E, 



Eo 



UJ2 



where use has been made of (overall) energy conservation, Ei 
a) + b) of Fig. 14 = (terms in iterated LS equation) 



E[ + E'2, and hence 
N 



(E2) 



(E3) 



cj| (Ei - Ep - E2 - Ui + ir]) 

The second term above in the static limit reduces to N/ {uJi tUg) ? which exactly cancels the 
contribution of diagrams c)-|-d). These exact cancellations persist also in the k)-l) as well as 
u)-v) type diagrams of Figs. H] and so that in computing their contributions we only take 
into account the subset of (twenty, see below) time-ordered diagrams of topology as shown 
in those figures. 

For the type j) contribution we find that the cancellation between irreducible and recoil- 
corrected reducible diagrams is only partial, and the result given in Eq. fl3.36p corresponds 
to taking into account only the irreducible diagrams illustrated in Fig. [15] (the same subset 
considered in the evaluation of type u)-v) above). However, the remaining irreducible and 
recoil-corrected reducible diagrams produce an additional contribution of the form 



e||iV,,(g) 



(t-1 X T2), (q X k2)j + [q x (cti x kg)] • 



0-2 ■ k2 



k2 
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/) 9) h) /) 




p) g) r) s) t) 



FIG. 15: Subset of time-ordered diagrams for the contribution illustrated by the single diagram j) 
in Fig. m See text for discussion. Notation as in Fig. [2l 



+ e^Mi) r-^^^ '^'^ "^ + 1^2. (E4) 

where the kernels Nij and Ni are 

which, however, does not lead to a Hermitian current density, since this would require 
j(ki,k2) = ki, — k2). We have ignored this contribution. 
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